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We investigate the parameter spaces of hybrid inflation model with special attention paid to 
the dynamics of waterfall field and curvature perturbations induced from its quantum fluctuations. 
Depending on the inflaton field value at the time of phase transition and the sharpness of the phase 
transition inflation can have multiple extended stages. We find that for models with mild phase 
transition the induced curvature perturbation from the waterfall field is too large to satisfy the 
COBE normalization. We investigate the model parameter space where the curvature perturbations 
from the waterfall quantum fluctuations vary between the results of standard hybrid inflation and 
the results obtained here. 



I. INTRODUCTION 

o 

' Thanks to vast data obtained from recent cosmological observations such as WMAP7 [1] inflation [2] has emerged 
as the leading theory of the early universe and structure formation. The simplest class of models of inflation is based 
'-"^ 1 on a scalar field, the inflaton field, coupled minimally to gravity [3, 4]. The inflaton field potential is flat enough 
to support long enough period of inflation to solve the horizon and flatness problems associated with the big bang 
cosmology. These simple models predict almost scale invariant, almost Gaussian and almost adiabatic perturbations 
which are very well consistent with the Cosmic Microwave Background (CMB) observations [1]. However, with the 
^3 ■ help of ever increasing accurate cosmic measurements different inflationary models can be discriminated based on 
i i ' their predictions for the power spectrum spectral index or the amount of gravitational waves. 

Recently there has been a revival of interest in hybrid inflation [5, 6], specifically the dynamics of waterfall field and 
whether or not significant large scale curvature perturbations can be produced during the waterfall phase transition [7- 
12], for earlier works see also [13-16]. In this work we search the parameter space of hybrid inflation in details and 
study the effects of the waterfall field dynamics carefully. We show that depending on the model parameters inflation 
can take place at different stages [17-19]. In conventional hybrid inflation [5, 6] inflation ends abruptly by the waterfall 
transition. Here we find situations in which inflation can proceed as in chaotic inflation [4] even after the waterfall 
transition. This is a new model of double inflation [20] in the context of hybrid inflation. Interestingly, as we shall see, 
in some situations inflation can actually have three extended stages. The key parameter is the ratio (f> c /Mp where <p c 
is the value of the inflaton field at the time of waterfall. With <j) c /Mp ^ 1 one can get a second period of inflation 
even after a sharp waterfall phase transition. 

We investigate the curvature perturbations in these models. The key constraint is that at the background the 
waterfall field ip is turned off so one has to introduce an effective trajectory \J (8ijj 2 ) from the quantum fluctuations 
of the waterfall field. This in turn plays crucial roles when one calculates the curvature perturbation power spec- 
trum. We exclude a large class of model parameters where the curvature perturbations cannot be normalized to the 
COBE/WMAP amplitude [21, 22]. 

The rest of paper is organized as follows. In section II we present our setup. In sections III, IV and V we consider 
a particular model in which inflation has three extended stages. After presenting the background solutions, we study 
the waterfall quantum fluctuations and the induced curvature perturbations in this model using the (5V-formalism. 
In section VI we investigate other limits of the model parameter space and calculate the curvature perturbations. 
We investigate the model parameters where the induced curvature perturbations from the waterfall float between 
the results obtained in sections III-V and the results of standard hybrid inflation. Some technical calculations of the 
waterfall field dynamics arc relegated into Appendices A, B and C. 
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II. HYBRID INFLATION AND ITS PARAMETERS SPACE 



As mentioned above we would like to investigate different limits of hybrid inflation potential parameter space and 
examine the effects of the waterfall transition on curvature perturbations. The potential has the form 

V-) = ^(V-^) 2 + + ^.9 W , (2.1) 

where 4> is the inflaton field, tjj is the waterfall field, and A and g are two dimensionless couplings. The potential has 
global minima at (0,-0) = (0,±M/\/A). 

In our analysis below, we retain some key properties of standard hybrid inflation as proposed in [5, 6] while relaxing 
other conditions. More specifically, as in conventional model of hybrid inflation, we keep the assumption of the 
vacuum domination. That is, during the first stage of inflation and before the phase transition starts, inflation is 
driven basically by the vacuum energy M 4 /4A. Secondly, we assume that the waterfall field is very heavy so it quickly 
relaxes to its instantaneous minimum ip — before the phase transition. Besides these two key assumptions, we relax 
the other conditions on model parameters. For example, in the original version of hybrid inflation 4> is sub-Planckian, 
while here we relax this condition and shall assume <j) c /Mp to be arbitrary where 4> c = M/g is the value of the inflaton 
field at the onset of waterfall. Furthermore, we keep the ratio g 2 /A arbitrary whereas in conventional models of hybrid 
inflation it is usually assumed that g 2 ~ A. Before we proceed we mention that similar studies have been performed 
in [23, 24]. These works mainly concentrated on the classical fields dynamics whereas in our analysis the quantum 
fluctuations of the waterfall field play key roles. 

The first period of inflation takes place at <f) > <f) c and ip = 0. As stressed above, as in a conventional model of 
hybrid inflation we assume this stage is vacuum-energy dominated, 

Mi 1 2 ,2 

» V , (2.2) 



with the Hubble parameter given by 



4A 



V((j),0) M 4 



H2 ~ 3Mj/ ~ 12XM 2 ' (2 ' 3) 

where Mp = l/y/8irG is the (reduced) Planck mass. We assume m 2 <C H 2 so that 4> is sufficiently slow-rolling during 
the first period of inflation before the waterfall. Also, as mentioned before, we also assume that the waterfall field is 
very heavy so that tp = is classically realized before the waterfall phase transition. These two assumptions about 
the masses of the inflaton and the waterfall fields can be translated into dimensionless parameter a and /? defined by 

m 2 M 2 

a ^Jp > (2-4) 

The fact that <j> is slow-rolling implies that a -C 1 whereas the assumption that ip is very heavy requires j3 ^> 1. 
The condition of vacuum domination, for d> near <t> c , can be written as 



f,2 
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Knowing that a<l this constraint indicates that there is no sub-Planckian restriction. Indeed one can also assume 
4> c >• Mp and a second stage of inflation may follow after a sharp phase transition. Alternatively, one can cast the 
vacuum domination condition in the form 

C=-^>1. (2.6) 
a A 

This condition indicates a hierarchy of the model parameters required to have the vacuum domination. However, the 
vacuum domination condition expressed by Eq. (2.5) shed more lights because it singles out the ratio (j) c /Mp as the 
key parameter which effectively distinguishes different limits of hybrid inflation as we shall see below. 

After <p < 4>c the waterfall becomes tachyonic triggering an instability in the system. We assume that this phase 
transition is mild enough such that it will take a few tens of e-folds for the waterfall field to settle down to its 
global minimum. Again this is possible if <j) stays near the critical point for a sufficiently long time. As we shall see 
below this condition requires m 2 M 2 <C H 4 . After the waterfall transition the quantum fluctuations of ip becomes 
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tachyonic and will grow exponentially. The accumulative effects of these growing modes, i.e. (Sip 2 ), through the 
interaction term changes the dynamics of the system. Soon after the gradient term in the waterfall held dynamics 
becomes negligible compared to the effective mass term the quantum fluctuations 5ip obey the same equation as 
the backgrou nd ho mogeneous ip equation. Hence they effectively behave like a homogeneous classical background, 
^classical = V '(Sip 2 ) , however see also [11]. 

The tachyonic growth of Sip continues until the self interaction term Xip 4 becomes important. We call this stage the 
waterfall stage. It ends when ip settles to its global minimum ip m in — ^j=- an d the symmetry breaking is completed. 
The other important effect of the tachyonic growth of dip is its back-reaction on the inflaton held via the interaction 
g 2 <p 2 ip 2 . This can increase the effective mass of the inflation held signihcantly ending its slow-roll sharply. The key 
question is which of the two back-reactions, Xip 4 or g 2 cp 2 ip 2 , becomes important sooner. In standard hybrid inflation 
with j3 extremely large, say /? <~ 10 3 , <p c <C Mp and g 2 <~ A these two back-reactions become important more or less 
at the same time. In this view the time of end of inflation is nearly the same as the time when the symmetry breaking 
is complete. 

In this work we distinguish the completion of symmetry breaking and the termination of inflation as two different 
phenomena. We will show that, depending on the model parameters, one can consider situation where the waterfall 
stage and symmetry breaking completion happens at an early stage of inflation, followed by a long period of the 
second stage of inflation. As we shall see the key parameter in the analysis is the ratio <p c /Mp. 

The effective mass of the inflaton can be read as 

m 2 eff =m 2 +g 2 (6ip 2 ) (2.7) 

Knowing that m 2 <C H 2 (or equivalently a « 1), inflation can only end when g 2 (5ip 2 ) ~ H 2 ~ M 4 /12XMp. Denoting 
this value of ip by iptrig which triggers the fast roll of <p towards its minimum, we find 



M 4> c 

^ tria ~ Yl\g 2 M 2 p ~ YlMy mm ' ! 2-S 1 



where ip m i n = For <p> c /Mp < \/l2, one has iptrig < ipmin so the back-reaction of ip on (p becomes important 

sooner than the self-interaction Aip 4 so inflaton held rolls quickly to its minimum before the symmetry breaking 
completion. During this transition, inflation may proceed briefly via the remaining vacuum energy but inflation will 
end quickly once the waterfall field nearly settles down to its minimum. This is the case in conventional hybrid 
inflation. On the other hand, for (p c /Mp > \/l2 one has iptrig > ipmin, meaning that the self interaction Xip 4 becomes 
more important than the back-reaction g 2 <p 2 ^p 2 . As a result, the waterfall field settles down to its minimum while 
inflation still proceeds afterwards as in chaotic inflation. Of course the dynamics of symmetry breaking is a violent 
phenomena which can cause a sudden violation of the slow-roll conditions. However, due to the attractor nature of 
the chaotic inflation potential, the slow-roll conditions are restored quickly for the second period of inflation. 

In our analysis it is convenient to use the number of e-folds dN = Hdt as a clock. For convenience we set N = 
at the time when the comoving scale corresponding to the present Hubble radius crossed the horizon during inflation. 
We denote the number of e-folds till the time when the waterfall field becomes unstable at <p = <p c by N — N c . Taking 
this epoch as the reference point we also define n = N — N c , so for the period before phase transition (after phase 
transition) one has n < 0(n > 0). 



A. 4> 2 C > V2M% 

For (pi > 12Mj>, one has iptrig > ipmin and the back-reaction of ip on itself becomes important first. The time 
when it takes for the waterfall field to settles to its global minimum depends on the sharpness of the phase transition. 
Denote by n — n p the number of e-folds it take for the waterfall field to settle down to its minimum. For a very sharp 
phase transition, n p < 1 whereas for a mild phase transition n p ^> 1. 

For n > n p , the effective potential for inflaton field is V e ff — \m 2 e ^(p 2 with m e // given by Eq. (2.7). For (p 2 > YIMp 
this leads to a period of chaotic inflation where (p evolves as 

<P 2 (n) = <p 2 p - 4M 2 (n - n p ) , (2.9) 

where <p p indicates the value of <p at the time of n p . As in chaotic inflation, inflation ends at <p = <pf where the slow-roll 
conditions are violated. The number of e-folds during the chaotic era, An c h 7 is 

An ch = %^ (2.10) 
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The total number of e- folds, therefore, is 

Ntot = N c + n p + An ch . (2.11) 

Here we have called the number of e-folds from N = until the end of inflation the total number of e-folds for 
convenience, but it should be noted that this number is nothing to do with the actual total number of e-folds of 
inflation which may well be almost infinite. 

To solve the flatness and the horizon problem we require N tot > 60. This picture only depends on the assumption 
that <j) 2 > 12Mp and is independent of the values of j3 and the ratio of the couplings g 2 /\. As we shall see below, 
the sharpness of the phase transition (the magnitude of n p ) depends on the parameter ~ a/3/3. If -C 1 then 
phase transition is mild and n p > 1. In this limit, as can be seen from Eq. (2.11), inflation has three stages. For this 
to happen, knowing that a <~ l/N tot , one requires that 1 <C (3 <C N tota i so (3 <~ 10 or so. On the other hand, if k$ > 1 
then the phase transition is very sharp and n p < 1. However, due to smallncss of a, to have k$ one requires that 
/3 » N tot , i.e. /? - 10 2 - 10 3 . 



B. <t> 2 c < Y1M% 

In this limit, i\)\ rig < ^ 2 min so the back-reaction g 2 ip 2 (j) 2 becomes important sooner than the self-interaction Ai/> 4 , 
triggering the inflaton field to roll down quickly to its minimum = 0. Thus inflation cannot continue for an extended 
period in the form of chaotic inflation, and one expects that the completion of symmetry breaking and the end of 
inflation occurs more or less at the same time. Specifically 

Ntot — N c + n tr ig + An^^trans + ^ n ip~trans ■ (2.12) 

Here ntrig denotes the time when the back-reaction of ip on the inflaton mass become significant, triggering its fast 
rolling to <j> = 0, whereas An^_t ran8 denotes the number of e-folds it takes for the inflaton field to move toward its 
minimum until the time when the coupling term g 2 (j> 2 ip 2 becomes negligible compared to the waterfall bare mass term 
—M 2 ip 2 , i.e. until when <j) 2 <C <f> 2 . Since <fi evolves like e -3 ™' 2 this period is rather short, An^-trans ~ 1. Once the 
value of <j) becomes negligible, the tp field becomes even more tachyonic and rolls quickly dow to its minimum like 

For large /3 this period of the transition (symmetry breaking completion) is very quick, An,/,_ trarls <C 1. One can 
show that ntrig ~ An^^trans ~ < 1, so inflation ends abruptly after the waterfall phase transition as in standard 
hybrid inflation. 

Having specified two major scenarios of hybrid inflation based on the ratio 4> C /Mp, we now proceed to the study 
of curvature perturbations in these scenarios. The case with <p 2 < 12MJ, is somewhat similar to standard hybrid 
inflation as studied in [7-11]. We shall come back to this limit in section VI. The case with cf) 2 > 12MJ, is the one 
that needs careful considerations. In the next three sections, as a sample example, we consider in details a model with 
cf) 2 > 12M|> and af3 C 1 so n p > 1 and inflation has three extended stages. This case perhaps is the most elaborate 
model in our case study which would be helpful in the classification of hybrid inflation models in section VI. Once 
we have performed the curvature perturbation analysis in this case, we are able to make a connection to the case 
(f> 2 < 12Mp and provide a global view of the issues of curvature perturbations and the waterfall dynamics in different 
limits of the parameter space of hybrid inflation. 



III. BACKGROUND DYNAMICS 



The analysis performed in this section and sections IV and V are for the case (f> 2 > 12M|> and a(3 <C 1. In this 
section we study the background dynamics in details which is also applicable to other cases. 
We assume the standard metric, 

ds 2 = -dt 2 +a(t) 2 dx 2 , (3.1) 

where a(t) is the scale factor. As mentioned before, we assume the vacuum energy dominance for most of the 
inflationary era. Hence we have 

a{t) cx e Ht , (3.2) 

where H is given by (2.3). 
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In this approximation, the dynamics of the system is solely described by the field equations for <p and ip. The field 
equations for the classical, homogeneous background fields are written as 

cp" + S<P' + (a + g 2 ^<f> = 0, (3.3) 

V + 3V + : 9 2 ^ + ^)^ = 0, (3.4) 

where the prime denotes the differentiation with respect to n, ' = d/dn, where n — N — N c , and the dimensionless 
parameters a and j3 are introduced in Eq. (2.4). We assume a<l and /? » 1. Here and in the following, we denote 
the classical fields by <p and ^ while the quantum fluctuations by Sep and Sip. 

At early times (— n) ^> 1, the classical waterfall field is trapped at ip = and the quantum fluctuations are 
completely negligible. Hence if we simply solve the classical field equation (3.4) with the initial condition ip = ip' = 0, 
ip would stay at the origin forever even at n > 0. However, because of the exponential growth of the quantum 
fluctuations of the waterfall field after the transition, the expectation value (Sip 2 ) becomes non-negligible and the rms 
value soon starts to behave as a classical field ip = \J (Sip 2 ). In fact we will show in the next subsection that (Sip 2 ) is 
dominated by the superhorizon modes, hence the rms of it is indeed observed as a classical background for an observer 
within each Hubble horizon region. In particular, in Eqs. (3.3) and (3.4), ip 2 in the interaction terms g 2 ip 2 /H 2 and 
Xip 2 /H 2 are to be considered as (Sip 2 ). 

Here we provide the solutions for the (p and ip evolution, starting with the dynamics of (p. Let us first assume that 
the interaction term g 2 ip 2 (p 2 is negligible. Namely, we consider the stage when a ^> g 2 ip 2 /H 2 in Eq. (3.3). With this 
assumption one can easily solve it to obtain 

<p(n) = </> c exp (— r n) (n < n p ) , (3.5) 

with 




(3.6) 

This stage ends when the value of g 2 ip 2 /H 2 becomes comparable to a. To simplify the dynamics, we assume that 
this happens at about the end of the waterfall stage, n = n p , when ip 2 = (Sip 2 ) levels off to the local minimum. 

Now we look into the dynamics of ip. Before the transition ip = but it effectively becomes non- vanishing after the 
transition. During the waterfall stage, we can neglect the self- interaction \ip 4 , and the equation for ip simplifies to 

ip" + 3ip' + fi(e~ 2rn - 1)V = 0. (3.7) 



The general solution is given by a linear combination of the Bessel functions J v (z) and Y v (z), where v = V / ^+ 9 / 4 anc ^ 

z = — ■ However, one can get a better view of the ip solution using the WKB approximation as we demonstrate 
below. 

Defining ip — vfe" 3 ™/ 2 , Eq. (3.7) is rewritten as 

(-^-2f3rn + 2f3r 2 n 2 ^j * = , (3.8) 

in which we expanded the exponential in Eq. (3.7) to second order in rn <C 1. Using the WKB approximation one 
has 



* <~ exp 



(l)i\/^ + 2/3rn-2Pr 2 n 2 



(3.9) 



As can be seen from Eq. (3.7) (see also Eq. (4.4) below) the sharpness of the phase transition is controlled by the 
product fir. The larger is fir the greater is the effective tachyonic mass of the waterfall during the phase transition, 
causing a sharp phase transition. This means that, by keeping the product fir small one basically reduces the effective 
tachyonic mass of the waterfall field during the phase transition so inflation can continue for an extended period 
during the course of the phase transition. In the analysis below, we are interested in a situation where fir <?C I so the 
phase transition is not sharp and n p 3> 1 . This is a new limit where the curvature perturbation shows some interesting 
behavior. However, for this to happen one requires /? < 1/r so /3 can not be arbitrarily large. For reasonable values 
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of r this requires j3 ~ 10. On the other hand, if one is interested in the limit where fir > 1, corresponding to n p ~ 1, 
then the analysis of the curvature perturbation from the waterfall is basically the same as in standard hybrid inflation. 
We will comment on this case later. In our analysis it is useful to define 



K * = T Si "9" <1 - 



(3.10) 



Calculating the integral in Eq. (3.9) to leading orders in and setting the initial condition at the moment of the 
phase transition n — yields 



il>{n) = -0(0) exp 



4 2 3 



(3.11) 



where 0(0) = \J (<50 2 (O)). The good agreement between this solution and the full numerical solution can be seen in 
Fig. 1. Note that the term in Eq. (3.11) containing fi^n 3 is much smaller than the term n^n 2 so we shall safely ignore 
it in most of our analysis below. 

The above solution is valid until the self interaction term Xtp 4 becomes important. Eventually the self-interaction 
becomes important and settles down to its local minimum because of the large induced mass, <9|V (0mm) = 
2(M 2 — g 2 (j> 2 ) ~ 2M 2 . Here note that the value of <j> will be a fraction of C by the time ip 2 reaches V'mm(^) ( see the 
next paragraph) and hence we may assume M 2 ^> g 2 4> 2 . Denoting the time at which ip reaches the local minimum 
by n — n p we have 



■0 2 (O)exp 



.3 „3 



M 2 

.w = - 



M 2 



(3.12) 



or to leading order 



2re„ 



■In 



A0 2 (O) 
M 2 



(3.13) 



This equation tells us how the total number of e- folds of the waterfall stage depends on tp 2 (0). As discussed in detail 
in the next section, it will play the central role in the evaluation of the curvature perturbation. 

Now, since we have assumed that g 2 ip 2 /H 2 becomes comparable to a at the end of the waterfall stage, we should 
examine its effect on the evolution of (f>. Inserting the solution (3.11) for tp into Eq. (3.3), one finds instead of Eq. (3.5), 



4>{n) = 4> c exp 



-rn 



gV(o) 

3H 2 



/ e 2K * n ' 2 dn' 
Jo 



(3.14) 



At the end of the waterfall stage, n = n p , we have ip 2 (n p ) = ip^in^) — M 2 /\. Hence the above gives 



4>(n p ) = <pc exp 



-rn„ 



3A 



.exp 



a 

-—fir, 



(3.15) 



where we have used Eq. (2.6) in the second equality. 

After tp 2 settles down to V'minW > ^- e - a * tne stage n > n p , we can plug this minimum value into the potential to 
obtain the effective potential for 0, 



V eff {<l>) = V(^ip mm (0)) = \m 2 4> 2 + ^M 2 <j> 2 - |^0 4 



m 2 + ^—M 2 
A 



(3.16) 



where the second line follows from our assumption that we have M 2 > g 2 (j) 2 by the end of the waterfall stage. Hence 
the g 2 (f) 2 /2 term may be ignored at the stage n > n p . The effective mass of </> at the stage n > n p is given by 



,eff 



m 



a 2 
X 



m 2 {l + C). 



(3.17) 



We note that the potential now behaves exactly the same as in standard chaotic inflation with m 2 — > m 2 (l + C). In 
other words, the universe is no longer vacuum-dominated at n > n p and inflation proceeds as in chaotic inflation. 
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n 



FIG. 1: Here we present the background fields dynamics in units of Mp from the phase transition n = to the end of inflation. 
The solid blue line and the solid red line, respectively, show the full numerical solutions of the <f> and ip evolutions with the initial 
value ip(n = 0) = (<5?/> 2 (0)). The dashed green line and the dashed brown line, respectively, show our approximate analytical 
solutions for the corresponding fields. The agreement between them is very good. In this picture, the time of phase transition 
completion is n p ~ 30. The numeric analysis are for the parameters, M = 5x 10~ 5 Mp, m — 2 x V0~ 6 Mp, g 2 = 2 x 10~ n and 
A = 2 x 1CP 9 corresponding to /3 ~ 8.5 and ~ 0.02. 



By using the above form of the effective potential, and assuming <p(n p ) ^> Mp, the equation for <p at n > n p reduces 



to 



-8M p Vn-n p ) = (0 2 -^)-(l + C- 1 ) 



In 



l-if(l + ^ 1 ) 



1 



which for 



bp < 4> c and C > 1 can be approximated to 

<t> 2 (n) =4$- 4M*t (n - n p ) (n > n p ) , 



(3.18) 



(3.19) 



where ^ p = (b(n p ) — 4> c e " lp . The slow-roll inflation ends when <f> becomes of order Mp as usual (approximately 
at (j) = 4>f ~ a/2/3 Afp). Here it should be stressed that the end of inflation is solely determined by the value 
of (j) alone. This is an important point to be kept in mind when calculating the curvature perturbation using the 
5 N- formalism [25] (see discussions in Sec. IV C). The overall evolution of the classical fields is depicted in Fig. 1. 

It is worth noting again that the above evolutionary picture holds at any horizon patch of the universe. Until the 
time of transition, i/j — J (Sip 2 ) cannot be regarded as classical. But as soon as the tachyonic instability sets in, it 
starts to behave like a classical field. Another important point is that the initial amplitude of these fluctuations at 
n = varies from one patch to another because of the very nature of the quantum fluctuations. As will be discussed in 
Sec. IV below, this implies the spatial variation of the number of e-folds from the time of transition to the end of the 
waterfall stage n = n p , which is fixed by the condition (8ip 2 (n p )) = ipmini'fi)' giving rise to the curvature perturbation 
at the end of inflation. It may be also noted that the quantum fluctuations of the waterfall field play the role of 
entropy perturbations if we regard <j> as the inflaton field. However, as we have described in the above, since they give 
rise to a non- vanishing classical background field ip whose evolution is crucial to the dynamics of inflation, the system 
is better regarded as a two-field inflation model after the waterfall transition. 
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IV. CURVATURE PERTURBATION FROM THE WATER-FALL FIELD 

In this section we study the quantum fluctuations of the waterfall field 5ip and calculate the curvature perturbation 
from them by using the JiV-formalism. 

A. Quantum fluctuations 

The equation governing fluctuations of the waterfall field in real space is 

dr + 3^' - ^V 2 ^ + (-13 + ,g 2 C + 3A H = , (4.1) 



where ip is the classical background. As we discussed in the previous section, we eventually identify it as ip 2 — (Sip 2 ). 
At early times before the transition or even after the transition, the self interaction term A(5ip 2 } / H 2 in the mass term 
of Eq. (4.2) is negligible until the end of the waterfall stage. In this approximation Eq. (4.1) in momentum space fe 
is given by 

HI + 3<5Vi + (J^p -f3 + g 2 ^y^ = 0, (4.2) 

where 

^-J^5iP(x)e- ik -* = (a k Mn) + al k Mn)) . (4.3) 

Here a& and a) k are the annihilation and creation operators, repsectivly, with respect to a suitably chosen vacuum 
and ■0fc( n ) is the positive frequency function. 

Substituting the solution (3.5) for <j>, one obtains 

H l + 3( ty£ + (^e~ 2n + /3 (e- 2rn -1)^5^ = 0, (4.4) 

where k c is the comoving wavenumber that crosses the horizon at the critical point n — 0, k c — Ha(n = 0). 

We assume a sufficiently small r such that r\n\ <C 1 for all scales k of cosmological interest. Hence we have 
/3|e~ 2r ™ — 1| ~ 2/3r|n|, which is always smaller than unity because of our assumption that k^, <C 1. Hence when the 
wavenumber k is inside the horizon, k > k c e n , one can neglect the term proportional to j3. The general solution in 
this limit is given in terms of the Hankcl functions H^) 2 (x) and H^ 2 (x) by 

S^ k = e- 3 "/ 2 (ciffg^x) + c 2 H^ 2 (x)) ; x = ^e~ n . (4.5) 

As usual, the natural choice of the vacuum is that it approaches the Minkowski mode function in the short distance 
limit, 

Sipk — > — y= as fcr — > — oo , (4.6) 
a\/2k 



where r is the conformal time, dr — dt/a. With this initial condition the positive frequency function for the modes 
inside the horizon is obtained to be 



W k (n) = J^e-^ 2 HV(±e-"). (4.7) 



4fc3 3 / 2 \k, 

At the time of horizon crossing n — n k when e Uk = k/k c , the above expression gives 



|<^(n fe )|=,-Jp. (4.8) 
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n 

FIG. 2: The amplitude of the mode function |-i/>fc(n)| as a function of the number of e-folds n is shown for a mode which leaves 
the horizon at n — 6. The solid blue line is obtained from the full numerical analysis whereas the dashed red line is from our 
analytical expression (4.9). All numerical parameters are the same as in Fig. 1. 



After horizon crossing one can neglect the spatial gradient term (k/k c ) 2 e~ 2n in Eq. (4.4) and the equation for tj} k 
becomes identical to that for the classical background, Eq. (3.7). Hence the solution takes the form (3.11). Matching 
it with (4.8) at n = n k , we obtain an approximate expression for ^ at n > n k as 



6ip k (n) 



H 



: exp 



«v(« - n k) - gM n ~ n k) 



(4.9) 



One can see a good agreement of this analytical expression with the numerical result in Fig. 2. 



B. Waterfall stage 



Given the solution to the mode functions, an important quantity to compute is the squared fluctuation Sip 2 , because 
it determines the duration of the waterfall stage. It is given by 



5ip 2 {n, x) = Sip (0, x) exp 
Sip 2 {0,x) -- 



d 3 kd 3 k' 



k.k'<aH (27r) 3 



<5Vfc(0)^l(0)e j(fe - fe,) - x . 



where 



^fc(O) = M0)A k = -JL= e - K *<+* K *<A k ■ A k = 



'2k* 



(4.10) 



(4.11) 



and we have absorbed an irrelevant phase factor in the definitions of a k and a_ k - Note that the integral is cut off at 
the horizon scale because the modes on subhorizon scales are just the standard vacuum fluctuations which are to be 
regularized to zero. 

The above expression for 8ip k shows that it is a classical random fluctuation, because it is proportional to A k which 
commutes with A^ p (= A_ p ) for any pair of k and p. Hence for each horizon-size patch one would observe Sip 2 as a 
homogeneous classical background which varies smoothly over scales larger than the horizon scale. Thus on a given, 
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sufficiently large scale, say the comoving scale of the present Hubble horizon size, one can calculate the mean value 
(5tp 2 (n)) and the fluctuation 

AiP 2 (n,x) = S^ 2 {n,x) - (S^ 2 (n)) , (4.12) 

where (Sip 2 (n)) determines the homogeneous background while Aip 2 (n, x) gives rise to the curvature perturbations 
on superhorizon scales. 

Let us first consider the background dynamics. Taking the expectation value of Sip 2 , we find 

(Syj 2 (n)) = W 2 (0)) exp[2/tyn 2 - jj K 2 n 3 ] , (4.13) 

where 

d 3 k 



w ' 2<0 » - L„ (0 lfc(o)|3 - (I) U t oxp| - 2k »"* + §** ■ 



(4.14) 



By changing the variable of integration to rik = ln(fc/fc c ), we may approximate the upper and lower limits of the 
integral to plus and minus infinity, respectively, for n^n 2 >• 1, and neglect the term n 2 ^ n\ in the exponent, to obtain 

-) j_Jn k *M-^nl) = -^—. (4.15) 

The waterfall stage ends when the self- interaction term becomes important and {Sip 2 ) reaches the local minimum 
i'mini^) given by Eq. (3.12). Thus we have 

W 2 K)> = W 2 (0)> cxp[2^n 2 - „3] = ^. n (0) . (4.16) 

Neglecting the k 2 n 3 term, we may approximately solve this for n p to find 



Before closing this subsection, let us make an important comment concerning the meaning of Eq. (4.11). If we 
would literally consider it as 5ip 2 (n, x) at n = 0, then it could only take account of the comoving wavenumbers k 
smaller than k c . However, as can be seen from Eq. (4.10), it takes account of all the comoving wavenumbers which 
has crossed the horizon at n — rife. This implies that 5ip 2 (0, x) introduced in Eq. (4.11) is defined just for convenience 
to express 8ip 2 (n 7 x) at the end of the waterfalls stage n = n p . 



C. (5iV-formalism 



We now evaluate the curvature perturbation from inflation by using the SN- formalism [25]. In this formalism, 
one first calculates the number of e- folds for the background homogeneous universe from a fixed final epoch t = tf 
backward in time to an arbitrary initial epoch t, N(t — > tf). The final epoch tf must be chosen such that at and after 
t = tf the evolutionary trajectory of the universe is unique. 

Then one considers the perturbation SN by identifying the final time slice to be a comoving (or uniform density) 
hypersurface and the initial slice to be a flat hypersurface. Naturally there will be fluctuations in the matter fields on 
the initial flat slice which produce SN. Then the resulting SN is equal to the comoving curvature perturbation on the 
final hypersurface 7Z c (x,tf) = SN(x). In the present model, we assume that the evolution of the universe is unique, 
i.e. there is no isocurvature perturbation, at and after the end of inflation. As before we use n as the time instead of 
t in what follows. Hence in particular, N = n— rif where n/ is the number of e- folds from the critical point until the 
end of inflation. 

There are a couple of points to be kept in mind when evaluating SN in the present model. The first point is about 
the dependence of N on the fields 4> and ip. Since we have to relate the fluctuations in these fields to SN, we have to 
express N in terms of cj> and ^. Namely, we need an expression like N = N((/)(n),')p(n)). 

Going backward in time from the end of inflation up to the end of the waterfall stage n p < n < nf, since the value 
of the waterfall field is determined by the value of the inflation, -0 2 = V'min^): w hh the identification ?A = \J (S?P 2 ), 
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N is a function of only <p, N — N(<p). Since the scales that leave the horizon after the waterfall stage n > n p are 
presumably too small to be of cosmological interest, we focus on the stage n < n p in the following discussion. 

Here the very important point is that the time when the waterfall stage ends, n = n p , depends heavily on the initial 
value of tp at the beginning of the waterfall stage, tp 2 (n = 0), hence on the amplitude of the fluctuations (Sip 2 (0)). 
This implies that N is a function of both <j> and ip 2 through its dependence on n p , 

N = N^<j>(n), ip 2 (n)j at n < n p . (4.18) 

The second point is about the choice of the initial hypersurface. In the present model the comoving curvature 
perturbation due to Sep is conserved on superhorizon scales. Therefore, instead of choosing the time of horizon 
crossing for each comoving scale, we can take the initial slice to be that at the time of transition n = but with 5(p 
given by that on the flat hypersurface at its horizon crossing as given by Eq. (4.21) below. Therefore for ^-dependence 
in N, we have 

S+N = N{<j> c + 5<t> c , VO - N((f> c , =N (t> (<j> + 5ct>^ <t> c ) - N+ty -> <j> c ) , (4.19) 

where N${<j> — > <p c ) is the number of e- folds from the time when the infiaton had a value cp until the time of transition 
when (f> = <p c , and Sep and S(p c are, respectively, the values of the field fluctuation at horizon crossing and at the time 
of transition. The contribution 5<pN may be approximated to linear order in 8<j> as it is the same as the conventional 
slow-roll single-field contribution which contains negligible non-Gaussianity. Noting 

——50 = — -5(f> = — 77 = -—, (4.20) 

d<p r (j> r 4>' r (f> y ' 

where we have used Eq. (3.5), we can evaluate S^N as 

1 6<j) c (x) 1 5<f>(n,x) e rn 5<j)(n,x) 
r <p c r <p(n) r cb c 

where we have used the fact that 5<j) on superhorizon scales has the same time dependence as the background <f>. When 
evaluating the spectrum, as usual, for a given scale k we choose n to be the time of horizon crossing n = ln(fc/fc c ) at 
which we have {5(g) = H 2 /{2tt) 2 . 

The contribution from Sip can be taken into account by taking the difference between the number of e-folds with 
5ip 2 and ip 2 = {8iP 2 ), 

S^N = N(<p c , 8^ 2 ) - N(<p c , (t/> 2 )) , (4.22) 

where Sip 2 and (Sip 2 ) are to be evaluated at the time of transition n = 0. Rigorously speaking, for scales with k > k c , 
it should be evaluated at its horizon crossing n = nk(> 0). However, as noted at the end of the previous subsection, 
we have conveniently included the modes k > k c also on the hypersurface at n — in 5ip 2 {0, x) by formally extending 
their evolutionary behavior back to n = as if they were already outside the horizon. Hence one may calculate 
5^,N(x) by just considering the fluctuations in 5ip 2 (0,x) without separating the modes into the two ranges k < k c 
and k > k c . A justification of this point is discussed in more detail in Appendix A. 

Equation (4.17) shows how the number of e-folds from the time of transition until the end of the waterfall stage 
depends on the magnitude of the initial fluctuations (<5V' 2 (0)). This implies that the total number of e-folds of the 
waterfall stage fluctuates from one horizon patch to another due to the fluctuations in 5ip 2 (0,x), giving rise to the 
curvature perturbation. Namely, we have 

1+ n p ) 1 2^n 2 p {(6tP 2 (Q)) J ' ^ 

where 5n p {x) is the spatial fluctuation in the number of e-folds from the time of transition until the end of the 
waterfall stage. Apparently we have S^N = 5n p (x). 

In the limit 8n p /n p -C 1, we can linearize the above to obtain 



(4.25) 



AK,^n p V W 2 (0)) 

One may further expand the above logarithm, by using ln(l + x) = x — x 2 /2 + • • • , to find 



5^N=—^— 

4:K^n p 



AiP 2 (0,x) 1 ^A^ 2 (0,a;)\ 



(8iP 2 {0)) 2 V (Sip 2 (0)} 
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where, following Eq. (4.12), we have Axp 2 (0,x) = Sip 2 (0,x) - (<5?/> 2 (0)). 

Here a comment is in order. One might worry that there could be a contribution to 8n p from the fluctuations in 
<j> through the ^-dependence of Vmin • Indeed there exists some contribution because the value of <fi at the end of the 
waterfall stage depends on the number of e-folds from the critical point as (j)(n p ) ~ <f) c e~ rnp . However, one can easily 
show that this contribution adds a term —r5n p to the terms inside the square brackets in Eq. (4.25), hence gives a 
negligible correction for K^n p /r <~ j3n p 1. 

After the completion of the waterfall stage the trajectory in the two-dimensional field space (4>, ip) is essentially 
unique, with ip being completely subject to </>, as can be seen from Fig. 3, and the end of inflation is synchronized on 
the uniform <fi hypersurface <fi — 4>f as in the case of single-field slow- roll inflation. This implies that there will be no 
contribution to the total SN from the era after the waterfall stage. 

As we mentioned before, the field trajectory is unique after the waterfall stage until the end of inflation. Setting 
n = nj and <p 2 {rif) — <p 2 in Eq. (3.19), the number of e-folds for this period An c h = rif — n p can be obtained as 

Hence the fluctuation in the number of e-folds from this slow-roll stage, S(An c f l ) is given by 

* (Anc ^§§ ^-"Sj^- (4 - 27) 

As is clear from this expression, 5{An c h) is suppressed by the factor r (<C 1) relative to Sn p . 

To summarize, the contribution to SN after the transition is equal to S^N to a good approximation. Hence the 
total SN is expressed as 

5N = 84N + 5i,N , (4.28) 

where 6$N is given by Eq. (4.21) and S^N is given approximately by Eq. (4.25). We note that, if we express Eq. (4.25) 
in the form 

S^N = N^A^P 2 + i^ 2>2 (A^ 2 ) 2 + • • • ,, (4.29) 
we have the relation among the coefficients, 

N "-^ N "'+-4^W>- (4 - 30) 
This is a useful relation when one wants to evaluate the non-Gaussianity in this model. 

V. POWER SPECTRUM AND CORRELATION FUNCTIONS 

This section is devoted to the calculations of the power spectrum and correlation functions of the curvature per- 
turbation. We find that the curvature perturbation induced from the waterfall is too large to be normalized to the 
COBE/WMAP normalization. 

A. Power Spectrum and spectral tilt 

The most important result of the SN formalism is that one can identify the curvature perturbation on the final 
comoving surface 7Z c (x), which is conserved thereafter, by the fluctuation in the number of e-folds from an initial flat 
surface to the final comoving surface, lZ c (x) = SN(x). Taking the Fourier component of Eq. (4.28) with the help of 
Eqs. (4.21) and (4.25), to leading order we obtain 

v = 1 A^1(0) 1 fkV SMnk) . , 

fe 4 K ^n p (<^ 2 (0)> + r \k c ) 4> c ' 1 ' ' 
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FIG. 3: Larger figure: Contour plot of constant N surfaces in the field-space are shown where ip and <j> are in Mp units. Solid 
blue lines correspond to different patches of the Universe with different initial amplitude of the quantum fluctuation 8ip 2 (0, x) 
at the critical point. The solid black curve denotes the trajectory with 8ip 2 (0, x) — {Sip 2 ). This figure shows the time interval 
between the critical point and few e-folds after the completion of the transition. Red contours show the constant TV contours. 
Counting the number of e-folds from the end of inflation, once the quantum back-reactions of Sip become important the constant 
N contours incline towards the constant ip lines while during the earlier stage of inflation the constant N lines coincide with 
constant cj> lines. Smaller (inner) figure: Trajectories of different patches of the Universe with different initial amplitudes of 
quantum fluctuations Sip 2 (0, x) at the critical point. This figure provides an extended range for the interval between the critical 
point and the end of inflation. The numerical parameters are the same as in Fig. 1. 



Using the fact that (A-0 2 ) fc = (<5^ 2 ) fc , the two point correlation function of the curvature perturbations can be read 



as 



<(<JV 2 ) fc (W 2 ) k ,) 1 (k\ 2r (6fa6<p k ,) n 



(ftfcft*) = T7> 2 2 77 , 2\2 + ~ F L . ( 5 - 2 ) 

where the second term is the contribution of the original inflaton field and the first term is the induced curvature 
perturbation from the waterfall field. 

To calculate the power spectrum of the curvature perturbation we first start off by calculating the correlation 
function of Sip 2 at the critical point, n = 0. As shown in [10], this correlation takes the form, 

(5ip 2 ) k (Sip 2 ) k ,) o = P sr {k){2^)H\k + k') 

d 3 q\S%(0)\ 2 Nfe-g|(0)| 2 S 3 (k + k') . (5.3) 



Here we estimate this integral which would suffice for our discussion here. In Appendix B we provide a more rigorous 
treatment of calculating this integral. Let us first estimate the contribution of the infrared modes in this integral, say 
\q\ < |fe|. The contribution of the long wavelength modes to the two point correlation function approximately can be 
found to be 



/' J d* q \5ip q (0r- (5-4) 

M<|k| 



The above integral counts the effects of large modes in the convolution. Note that we take the term |i5-0|fe_q|(O)| 2 out 
of the integral which is possible for small q. To simplify the analysis, one can take the upper bound to be k c so the 
contribution of the infra red momenta in the convolution can be read as 

P$ = \HM\ 2 |J f dn q e- 2 «*< = i(^ 2 (0))|%(0)| 2 . (5.5) 
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Now we estimate the contribution of the small modes in the convolution of integral in Eq. (5.3). By the small 
modes we mean the modes for which, |q| > |k|. By this assumption, one can simply find 



pUV _ 
^8xp 2 — 



1 



(2tt) 3 



/ 



d 3 q\SM0)\ 4 



H 4 



167T 2 fc3 



\g\>W 



In k/k c 



dn q e 



-3n q —An^,n q 



(5.6) 



As we are interested in the |<7| > |fc| to lower bound simply taken to be k = q in r.h.s. The second order term in the 
exponent suggests a natural cut-off n q < 1/^/2/c^, for the integral so one can simply ignore the quadratic term of n q 
in the exponent to obtain 



pUV 



<^ 2 (0)>|^ fe (0)| 2 . 



(5.7) 



One readily finds that the ultraviolet contribution from Eq. (5.7) is suppressed compared to the infrared contribution, 
Eq. (5.5), by a factor ^JK^, <C 1 so the power spectrum from the ip contribution mainly comes from the infrared modes 
and P s ^2 ~ P/^ 2 . 

Finally from Eq. (5.2) the curvature perturbation spectrum per unit logarithmic frequency interval, 



*-5?fl«-5Fl*.l" 



(5-8) 



is obtained to be 



" k 3 |^ fe (0)| 2 




2rn k 



2rn k 



4ir 2 r 2 /3 



(5.9) 



We need g <C 1 so the contribution from the infiaton field, V$ is sub-leading compared to VvJ R ■ 

With the curvature perturbation power spectrum dominated by Vn R the spectral index of the curvature pertur- 
bation can be easily found as 



1 



d\nk 



-4K^,rik = — 4K^,ln(fc/fc c 



(5.10) 



Since the curvature perturbation is nearly scale invariant but it has a slight running. It runs from blue for 

modes which exit the horizon before the transition to red for modes which exit the horizon after the transition. One 
can easily quantify this running as 



dn s 
din k 



-4k 



(5.11) 



This should be compared with the standard hybrid inflation where it is found that n s ~ 4 [7-11], a very blue-tilted 
spectrum. This is the key property which renders the curvature perturbations from the waterfall field harmless in 
standard hybrid inflation. 

Now we face the fatal problem in our model. We see from Eq. (5.9) that the curvature perturbation power spectrum 
cannot be normalized to the required COBE/WMAP normalization V-r, ~ 2 x 10~ 9 . With and ^K^n p <~ 1 one 

finds that P-r IR is many orders of magnitude bigger than the COBE amplitude. In comparison, in standard hybrid 
inflation this problem does not show up. As shown in recent analysis [7-11] the curvature perturbations from the 
waterfall field in standard hybrid inflation with (f> c <C Mp and k$ >• 1 scales like (k/k c ) 3 ~ e~ 3Nc which is totally 
negligible on the cosmological scales. 

One may think under what conditions the curvature perturbations induced from the waterfall can be kept sub- 
dominant on the cosmological scales. To address this let us first estimate the minimum e-foldings required between 
the transition time and the time of horizon crossing of the large scale (CMB) modes. We would like to impose the 
following constraint on the amplitude of these modes 



IR {k CMB ) < 2 x 10- 



(5.12) 
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in which k CMB is the momentum of cosmological sizes observed on CMB. By using Eq. (5.9), for typical values of 
our parameters, this yields 




\n k cMB | > 23 x , - — . (5.13) 



To have 50 — 60 e-folds after the transition, we must have n p + An c h = 50 — 60. But if An c h dominates it is essentially 
equivalent to a chaotic inflation model which is not interesting. Hence let us assume n p > 50 or so. Then Eq. (5.13) 
above implies that this is possible only for n^ MB < 0. But then the number of e-folds form the CMB scale till the 
end of inflation becomes too large. Hence it is impossible to hide the curvature perturbation from the waterfall field 
for the range of the parameters we consider unless the final chaotic inflation stage dominates. 

Now comes the natural question that what is the source of the problem? The key to this questions is the lack of 
a classical background ip trajectory so the quantum fluctuations Sipkin) determine the effective classical trajectory 
y/ (Stp 2 (n)). The ratio of quantum fluctuation 5ipk to this effective background trajectory represents the amplitude of 
the quantum fluctuations. This ratio in the coordinate space is 

Jity, (5.14) 



(2tt) 3 (<5V 2 (0)) 



which is not small enough to be normalized to the observations. To see this better compare the situation with the 
standard chaotic models. In models of chaotic inflation, the fluctuations are given by 5(j)/4> in which 5(f) <~ H while 
the value of the background (p field is another free parameter which can be tuned independently to get the right 
normalization. In other words, in chaotic models there is a well-defined background classical (p field which exists 
independent of the quantum fluctuations 5(p. However, in our model there is no classical background ip field. Instead 
the effective background waterfall field yj (Sip 2 (n)) is constructed from the quantum fluctuations 8ipk{n) and as a 
result there is no free parameter to tune the relative ratio k 3 \5ipk(fy\ 2 / (Sip 2 (0)) to a very small value. 

Having this said, one can naturally ask why this problem does not appear in standard hybrid inflation model. As 
mentioned above, the main feature in standard hybrid inflation is that the large scale wa terfall perturbations are 
highly blue-tilted, scaling like (k/k c ) 3 . On the other hand, the main contribution into \J (Sip 2 ) comes from small 
scales, modes which become tachyonic but remained sub-horizon till end of inflation. As a result the cosmological 
large scales, corresponding to the current Hubble size, are suppressed compared to these small scales by the ratio 
(k/k c ) 3 <~ e - 3Ar <=/ 2 and therefore are completely negligible in the curvature power spectrum. In contrast, in our 
model the power spectrum is nearly scale invariant, as can be seen from Eq. (5.10), so the large scale modes are not 
significantly suppressed compared to the small scales and both of these modes contribute equally into yj {Sip 2 ). As a 
result, one cannot tune the ratio k 3 \5tpk(0)\ 2 / (Sip 2 (0)) to a very small value. 

In passing we note that in hybrid inflation with a sharp phase transition if one reduces N c so the phase transition 
happens early during inflation, i.e. N c ~ 5, then the large scale modes may survive the suppression. We address this 
case in the next section. 



VI. GENERAL PARAMETER SPACE OF HYBRID INFLATION 

In last three sections we have focused on the model where (p 2 > 12Mp and Ky, < 1 so inflation has three extended 
stages where N to t = N c + n p + An c h- We have particularly studied this limit of the parameter space in great details 
for two reasons. First, this limit was not studied previously in the literature, namely the effects of the waterfall phase 
transition on curvature perturbations. Secondly, this limit of parameter space provides us a specific example in which 
one finds an unexpected result that the contribution of the waterfall field to the power spectrum is too big to satisfy 
the COBE normalization. As we have seen in the previous section, this is because the effective background trajectory 
(Sip 2 ) gets more or less equal contributions from both the large scale and small scale quantum fluctuations Sipk ancl 
as a result there is no free parameter to suppress the ratio \5ipk(0)\ 2 / (5tp 2 (0)) . 

For the other limit of parameter space, we have standard hybrid inflation with <p c < Mp and k$ > 1 so inflation 
ends quickly after the phase transition. As shown in [7-11] the induced curvature perturbation from the waterfall field 
is highly blue-tilted and is completely negligible on large scales. We would like to search the parameter space of hybrid 
inflation in general and examine the dynamics of waterfall phase transition and its contribution to power spectrum. 
We would like to see under what conditions the results vary between the results of standard hybrid inflation and the 
results we found in the previous sections. 

As mentioned in section II, two major classes of the parameter space were determined by the ratio (p 2 /12Mp. In 
the previous sections we have studied the case <p 2 > YIMp and Ky, < 1 in details. Now it is time to consider other 
limits of the parameter space. 
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A. 4> 2 C > YIM 2 P , > 1 

Unlike the case studied in the previous sections, the condition n^, ~> 1 results in a sharp phase transition similar 
to standard hybrid inflation. However, in contrast to standard hybrid inflation in which inflation ends shortly after 
the waterfall, in this case with fa 2 , > 12MJ, inflation continues for an extended period after the transition. On the 
other hand, knowing that a ~ l/N tot <C 1, to have » 1 one requires that /? ^> N to t, i.e. fi ~ 10 2 — 10 3 . As in 
previous sections, inflation has three stages and N tota i = N c + n p + An c j. Here we first briefly review the background 
dynamics and then compute the final curvature perturbations using the 5N formalism. 

Stage. 1 The first stage corresponds to time before phase transition, i.e. N < N c when the inflaton field reaches to the 
critical value, cj> = fa, after which the waterfall field becomes tachyonic. Fluctuations of the inflaton field, 8(f), 
which exit the horizon before the critical point cause variations on the duration of this period, SN C . By using 
Eq. (3.5), one has fa — 4>ie~ rrii , in which fa- and rij respectively are the value of the inflaton field and the time 
of horizon crossing for a specific mode. So the variation of iV c can be read as 



— ( -) (O.D 

fe<fe c 



r 



Since the next stage of inflation starts with <f> = fa, for each patch therefore the information associated to these 
perturbations is encoded in SN C 

Stage. 2 This stage is defined as the period after the transition until the completion of transition, n p , when tp becomes 

M 2 n 2 

W 2 )=^ min (<t>) = — -^<t> 2 . (6-2) 

As in standard hybrid inflation this period is very short, n p < 1, but since the inflaton field is super-Planckian 
inflation can resume for the third stage in the form of chaotic inflation. As shown in [8-11] the spectrum of the 
waterfall fluctuations which become tachyonic shortly after the transition is highly blue-tilted while for modes 
which become tachyonic later on the power spectrum is highly red-tilted. For modes in between the spectral 
index is running significantly from blue to red so one should avoid these modes in contributing to the total power 
spectrum. As we discussed before the amplitude of cosmological large scales is suppressed by the f actor e ~ 3Af ^/ 2 
compared to the small scales. Since the small scale modes have the dominant contributions into \J (Sip 2 ) as the 
effective classical trajectory, the contribution of the large scales is suppress in the power spectrum by the ratio 
e~ 3Nc . Below we will provide a bound on N c in which one can safely neglect the contribution of large scale 
curvature perturbations. 

Stage. 3 This stage is the period of chaotic inflation after the transition, n > n p . The end of this period is defined by an 
explicit value of the <j> field. The ip field in this period does not play any role because it is held on its minimum. 
The inflaton fluctuations which exit the horizon after critical point can change the number of e-folds during 
chaotic period. From Eq. (4.26) one has 

nf-n p - 4Af2 • (6.3) 
pi 

To calculate curvature perturbations using 6N formalism first concentrate on the fluctuations which exit the horizon 
before the critical point. One can borrow the results from Section IVC. The only difference is the dynamics of the 
waterfall field. In general suppose that the waterfall field has the following form after the transition 

^ 2 (n) = V> 2 (0)e 2/ K (6.4) 

For our case at hand, in which the phase transition happens similar to standard hybrid inflation [8, 10], we have 

/(«) - §V^7^ 3/2 ■ (6-5) 
The time at the end of transition can be found by using the equation, 

1& in = <ty 2 (0)e 2 '<"'>. (6.6) 
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As before, ignoring the slight dependence of 4> m in{4>) on <fi, the variations of the phase transition duration to the first 
order can be read as 



-1 AV> 2 (0,x) 

"^' ifw WW ' (6 - 7) 



in which AiP 2 (0,x) = 5iP 2 (0,x) - (8^(0)) . 

Adding up the contributions from the inflaton field, Eq. (6.1), and the waterfall field, Eq. (6.7), one finds the final 
curvature perturbation for modes k < k c , 

Finding the curvature perturbation for modes which leave the horizon after waterfall, k > k Cl requires some careful 
considerations. We relegate the details into Appendix C and quote the final result, 



2f(n p ) (srm 



f'(n P ) 



k>k c . (6.9) 



As is seen from this result, the curvature perturbation spectrum is smoothly taken over from 5(j> to Sip 2 as n increases 
from to rip. 

As mentioned before, during the second stage the power spectrum can run significantly from blue to red. To avoid 
this one has to make sure the phase transition happens late enough during inflation so the contributions of these 
modes are negligible on large scale power spectrum. Following the analysis in [10] (for a method different than SN 
see also [8]) the contribution of the waterfall field in total curvature perturbation is 1 

p K W),^ M . (fU0) 



i^p V k r 



The cosmological large scales relevant for CMB exit the horizon about 60 — 53 e-folds before the end of inflation, so 
one should check that P n ( s ^ is smaller than the observed COBE normalization Vn <~ 10~ 10 . This yields 



which leads to the following constraint 



K - 11 /3 e -3(JV e -7) <10 -10 > (6U) 



N c > 14.7- H l nKv ,. (6.12) 



The above relation indicates that if one takes n$ large enough then Vtz^^ is so suppressed that the transition can 
take place for N C >1. Below we study each stage in brief. 

B. eft < V2M% 

As discussed in section II inflation has four stages N tot = N c + n trig + An0_ trans + Arty —trans in which the <j) 
transition period is very short. 

Stage. 1 This stage, corresponding to N < N c , is completely similar to the case studied before. 

Stage. 2 This stage corresponds to the case where the back-reaction of ip on <f> becomes important triggering it to move 
towards its global minimum = 0. The value of the ip at the time of ntrig has no dependence on inflaton field 
so fluctuations of the inflaton field does not affect the duration of this period. Furthermore, one should now 
replace n p in previous subsection by n tr i g because the transition of ip shuts off at n tr i g instead of n p . One can 
show that 

2 1 i I Ymin 4>c \ ,„ -. \ 



1 Note that the notation a used in [10] can be translated to our kj, here as a = ^/6k^ 
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As in the previous sub-section suppose after the phase transition ip ~ Tp c e^ n ^ where the form of f(n) depends 
on the sharpness of the phase transition. Following the same steps as in Section IV C, the contributions of ip 
fluctuations into the variation of the waterfall transition number of e-folds is 

Sn ^ = WM ln {jmo)))- (6A4) 

Regardless of the details of the transition, one deduces that the contribution of the waterfall field quantum 
fluctuations into the total curvature perturbations is 

W = 0) = S „„, = ^||2|. (6.15) 

Stage. 3 This stage corresponds to the short interval An^- tr ans when (f> rolls quickly to its minimum. One can still 
consider the vacuum domination in this period as the waterfall field is far from its global minimum. The 
equation of the inflaton field can be read as follows 

(/)" + 30' + uj 2 {n)(j) = (6.16) 

in which 

w 2 (n)=a+ ^2 (fu7) 

and ip ~ iptrig^^- Defining <f> = e~ 3 / 2 ™$ one has 

^.+ u 2 (n)j $ = 0. (6.18) 

Shortly after the triggering time the term containing ui 2 becomes large and completely dominates over —9/4 
so one can use the WKB approximation to find the time dependence of the inflaton field during the fast roll 
towards = 0, 

<f>~ e -3n/2 e ±i/" W ( n ')dn'_ (g 19 ) 

This indicates (p scales like e~ 3n / 2 during its fast roll towards its minimum. The transition of 0, triggered due to 
back-reaction of the waterfall field, can be considered to be complete when the interaction term g 2 (p 2 ip 2 becomes 
negligible compared to the bare mass of the waterfall field 



a2 



« 4> 2 C . (6.20) 



Since the inflaton field slow rolls before the triggering point, its amplitude there is about <f> c . One can readily 
concludes that the fast roll of inflaton field completes in 

An^-trans ~ 1 • (6.21) 

So one can consider that this transition is somewhat instantaneous. 

Stage. 4 Soon after settles to its global minimum — the waterfall field experiences the large tachyonic mass — M 2 . 
One can check that in this period ip starts growing like ip(n) ~ e ^ n towards its global minimum. 

Again, to find the total curvature perturbation one should add up different contributions to 5N from different periods 
specified above. If 3> 1, corresponding to a sharp phaser transition, the situation is the same as in standard hybrid 
inflation. Since the inflaton field is sub-Planckian then inflation can not continue after the waterfall in the form of 
chaotic inflation so one basically requires N c ~ N t ~ 60. Furthermore, the large scale curvature perturbations do not 
get contribution from the waterfall field fluctuations which are suppressed by (-^-) 3 ~ e~ 3Na <C 1. Note that with the 

assumption > 1 we have /3 » N to tai i- c - P ~ 10 2 — 10 3 . 

However, if we assume « 1 so the phase transition is mild, then we face the same problem as discussed in 
section V, that is the curvature perturbations cannot be normalized to the COBE normalization. This can be seen 
from Eqs. (6.14) which has the same structure as Eqs. (4.24) and (4.25). 
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VII. SUMMARY AND DISCUSSIONS 

In this paper we have looked into the parameter space of hybrid inflation. We have retained two key properties of 
the standard hybrid inflation scenario, the vacuum domination condition and the condition that the waterfall field is 
very heavy during early stage of inflation. The latter condition is equivalent to /3 ^> 1 causing the waterfall field to 
quickly rolls to its instantaneous minimum ip = 0. Besides these two conditions, we have relaxed other conditions on 
the hybrid inflation parameters and searched for the limits of the parameter space where one can obtain curvature 
perturbations with the right amplitude. As we have seen, two key parameters in our classification are <p c /Mp and 
k$ where the latter parametrizes the sharpness of the phase transition. For a sharp phase transition, as in standard 
hybrid inflation, one has k$ 3> 1, whereas for a mild phase transition /j^<l. 

We have shown that models with (p 2 > 12Mj> and < 1 do not work because they can not satisfy the observed 
COBE normalization. Note that to have « 1 while keeping the waterfall field heavy requires that [3 ~ 10 or so. In 
this model inflation has three extended stages, the first stage is before the waterfall, the second stage starts when the 
waterfall rolls mildly to its minimum while inflaton is slowly rolling. The final stage is similar to conventional chaotic 
inflation models. As we have seen, the dominant contributions into the curvature perturbation power spectrum come 
fr om the w aterfall fluctuations. On the other hand, there is no background classical ip trajectory so one identifies 
\J (Sip 2 (n)) as the effective background trajectory. As a result, one finds that there is not a very small free parameter 
to tune the ratio (fc 3 / (27r) 3 ) 1(5-0/= (0) | 2 / (<5'i/> 2 (0)) ~ ,Jk^ to the small observed value. This is in contrast with, say the 
m 2 2 chaotic model, in which 5<p ~ H while <f> is kept as an initial value free parameter so one can tune 5 </>/</) into an 
arbitrary small number. On the other hand models with cf) 2 c > YIMp and k$ » 1 can be allowed. The contributions of 
waterfall field fluctuations are the same as in standard hybrid inflation. However, one has to satisfy inequality (6.12) 
in order to avoid a significant running of spectrum from blue to red on the CMB scales. 

Models with (p 2 < \2Mp and k$ ^> 1 are similar to standard hybrid inflation and one can satisfy the COBE/WMAP 
normalization. In these models inflation end abruptly after the waterfall phase transition. For having ~> 1, one 
requires that (i <~ 10 2 — 10 3 . However, in models with <p 2 < 12MJ, and k^«1 one encounters the same problem as 
above that the COBE/WMAP normalization cannot be satisfied. This corresponds to /? ~ 10 or so. 

After these discussions, one concludes that models of hybrid inflation with the assumptions of vacuum domination 
and a mild phase transition do not work because one can not normalize the curvature perturbation power spectrum 
to the required COBE/WMAP normalization. This corresponds to models with k^«1 and (3 <~ 10 or so. 
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Appendix A: Justifying integrating out sub- horizon modes in Eq. (4.10) 

Here we justify our approximation in integrating the amplitudes of all of quantum modes regardless of being super- 
horizon or sub-horizon at critical point, given by Eq. (4.10). In order to employ the 5N formalism we have to start 
with a well-defined background classical trajectory [26]. In our model, there is no well-defined background trajectory 
for the ip field. However, the accumulative effects of quantum modes Sip defines the effective ip background as \J (dip 2 ). 
To find this average we should consider the modes which are at any time super-horizon. In other words, we renormalize 
the two point function with cut-off defined by the comoving momentum A = k/a = H. By the above discussion the 
background ipb at any time can be read as follows 

= {8 ip 2 ) ren . (n) = jf ^ | e 2 ^" 2 . (Al) 

By using Eq. (4.9) one can simplify the above integral 

tt2 rn 

ip 2 b = ^ J dn k e- 2K ^ e 2 ^" 2 (A2) 
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In the above integral, at any given time time, only the super-horizon modes are considered. The above integral can 
be find in terms of error-functions 



(A3) 



By using the properties of the error-functions, one can verify that for the time n > 1/ y/2n^ the above expression 
approaches quickly to 



<**"<«>» 



(A4) 



Now let us compare with n p : the bigger is the value of ^/K^n p > 1 the better gets our approximation Eq. 

(A4). Let us now have an estimate for ^fK^n p . From Eq. (4.17) we have 



K^n 2 



In I -./•-!')-■! 



^32^ 



(A5) 



The term inside the square root is basically at the order of unity, therefore n^n 2 ~ — In >/X. For A sufficiently 
small one can arrange for K^n 2 3> 1. For example in our numerical analysis, A ~ 10 -9 which gives n^n 2 ~ 21 so 
our approximation Eq. (A4) is very good. Therefore for times well after the phase transition but before the phase 
transition completion one can use the approximation (A4) as a good estimate for the background ip field. 

In the SN formalism the evolution of curvature perturbations is considered only for super-horizon modes [25, 26]. 
Working on super-horizon scales means smoothing out the sub-horizon perturbations in a patch. At any time super- 
horizon modes are the modes with the comoving wavelength larger that the Hubble radius, k < aH. Therefore one 
can integrate the curvature perturbations by considering only the super-horizon modes at a given time. For example 
in our model, at the time of phase transition there are many modes which are sub-horizon, but as time goes by 
these modes become super-horizon and should be considered in the background evolution. This point is what we 
considered in evaluating the average trajectory oiip held. Therefore to calculate 5tp 2 (x), defined on every patch, only 
the super-horizon modes should be taken into account. This qualitative discussion verifies that for each patch and for 
times well after the critical point one can use the approximation 



5ip 2 (x,n) ~<5?/> 2 (x,0) e 2K <*" , 




(A6) 



The above arguments support our discussions at the main text about the main n-dependence of Sip 2 for every patch. 



Appendix B: Careful treatments of the waterfall induced curvature perturbations 



Here we study in details the integral in Eq. (5.3) to find the induced curvature perturbations from the waterfall 
field. 

From Eq. (4.9) for the mode function 5ip p , we have 



l<ty|fc-,|(0)| s 



Ml 

2fc? 



exp 



-3n| fc _q| - 2K^nf k _ ql 



(Bl) 



where 



r k 2 + q 2 - 2qk cos ( 



k 2 



(B2) 



and = (3r/3. 

The momentum integral contain an angular part plus radial part. Taking account of the azimuthal symmetry, we 
can write 



d 3 q = / 2irq 2 dq / ci(cos 9) 



(B3) 
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where cos 6* = k ■ q/kq. Changing the integration variable from cos 6* to n\k-q\, the angular integral becomes 

fj2 r nk +i 



d(cos0)|<ty| fc -,|(O)|' 



2k c qk J n]k _ q] 

The result of the above integral is in the form of the error functions 



/•nk+q 

/ dn exp [— n — 2/^n 2 ] . 

Tl It- — n\ 



/ 



d(=os S )|^ |t _„(0)|» = £-J^e 



,1/(8**) 



Erf 



1 + ^K^Uk + q \ 

2^2^ J 



Erf 



1 + 4K^n\ k _ q \ 



where 



and the error function is defined by 



n k+q = In 
n\k- q \ = In 

Erf(z) = 



k + q 

k c 
\k-q\ 

kr 



Jo 



'dt. 



(B4) 

(B5) 

(B6) 
(B7) 

(B8) 



Plugging Eq. (B5) into the total momentum integral one has 

<W' 3 ) t («</■% > = * 3 <* +<=')* 



Erf 



1 + Anipn\ k _ q \ 



(B9) 



The error function, Erf(z) varies from — 1 for z — > — oo to +1 for z — » oo. This transition is sharp and take place 
near z — 0. One can see that the overall behavior of this function is similar to tanh(z). For k not extremely small, 
the argument of the first error function in the square brackets is always much greater than unity, hence this term is 
almost 1. The second term depends on the difference between the amplitudes of two momenta, \k — q\, and can change 
from — 1 for q = k to 1 for q completely different from k. There is two criteria, first, a very narrow region in which 
the argument of the second error function become small or negative, in this region the difference between two error 
functions is considerable but the domain of the integration is very small. Second possibility is when the amplitude of 
the integrand is not too large but the domain of integration is large. This point is more clarified in the following 

We showed in the draft that the contribution of the modes in the first criteria is suppressed by a factor e^ 1 / 8 ^ as 
follows on the first contribution to the integral. 



2(2tt) S (k + k' 



IP 

4tt 2 V 2k^ 



(BIO) 



It's clear from the steps of the calculation that this suppression is due to narrow width of the integration domain. 
It's time to calculate the second contribution. This contribution comes from the region in which the arguments of two 
error functions in the big bracket are not too far apart and hence the difference between the error functions is small. 
As we pointed out before, the argument of the first error function is large and it's amplitude of is approximately equal 
to unity, outside the first regime discussed above the amplitude of the second error function is also large so we can 
use the following estimate for the behavior of the error functions with large argument 



(1 - Erfz) 



(Bll) 



Let us first clarify a tiny point. We defined before the first contribution to the integral comes from the region in which 

1 + 4/tyTl|fc_ g | 



< z bw 



and here we are going to concentrate on the region in which 



1 



> 1 



(B12) 



(B13) 
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one can say as Zb w — 1/4 you may miss some of contribution to the integral but we claim that this error is not 
considerable. By using the Eq. (Bll) one has 

Erf { 2^ J - Erf [ 2^ ) * — l + 4^n, fc _/ XP H*-«l - 2 ^ n l*-.l) 

One can also more simplify the above equation by ignoring the term +4K^n\k- q \ in the denominator of the fist fraction. 

This approximation can be verified by noting that the total term has as exp (%Kil>n 2 k _ q ^ dependence. 

By using the above approximation one has for the second contribution to the integral of the two point correlation 
function 

((# 2 ) fc (6^) k ,) 2 = S 3 (k + k')x4nH 2 J g|^(0)| 2 exp(-n |fe _ 9| -2^nf /t _ 9| ). (B14) 
\k — q\ 

By noting that n\ k _ q \ = In — , one has 

dq = ±k c e n \ k -^dn\ k _ q \ (B15) 

respectively for q > k and q < k. Again we claim that as this combination of the error functions are highly peaked 
around the q = k we take the explicit g-dependences out of the integral and identify q in them by k. By this trick 
and noting that the remaining integrand is even one has 

< (^ 2 ) fc (S^ 2 ) k , ) 2 ~ 5 3 (k + k') x SttU 2 |^ fe (0) | 2 J dn ]k _ ql exp (-2/^nf fc _,,) . (B16) 

By performing the integration finally we have 

( (^ 2 ) fe (<^ 2 ) fe/ ) 2 ~ (27rf6 3 (k + k') x i(^ 2 (0))|^(0)| 2 (B17) 

This relation is in complete agreement with our previous estimate for the two point correlation function, Eq. (5.5) 



Appendix C: Curvature perturbations for k > k c with (f 2 c > YIM^k^ >• 1 

In this appendix we are going to study the contribution of the quantum fluctuations with k > k c to the final 
curvature perturbations for the case 4\. > 12M 2 , using the SN formalism. For this purpose, we trace back the number 
of e-folds, N, from the end of inflation until time of horizon crossing of some specific mode during transition n. To 
avoid confusion, in this appendix we denote the number of e-folds counted backward from end of inflation by N. 
Namely, N = n — rif where n/ is the number of e-folds from the critical point <f> = 4>c until the end of inflation. Our 
strategy is to express Af in terms of the fields (f>(n) and ip 2 {n) (smoothed on every Hubble patch). 

As discussed in sections II and III, in the models with cf) 2 > 12M 2 , the transition is followed by a chaotic period of 
inflation. So the end of inflation is determined by the value of <f> alone as 

4> = <j> f w y/2/3M P . (CI) 
From this up to the end of the waterfall transition N is given by 

4> 2 - 4> 2 f d> 2 

"=^-3^ M - M '- (C2) 

where Af p is the value of Af at the end of waterfall transition. As shown in Eq. (3.12) the value of ^ 2 at Af = Af p is 

,2, s ,2 , , M 2 g 2 m P ? M 2 - Ag 2 N p M 2 

4> (n P ) = Tpmin(n P ) = 1 = 3; • (C3) 

Here Af p is given in terms of the number of e-folds from the critical epoch <f> = <f> c to the end of waterfall transition, 
which is denoted by n p , as 
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Hence ip 2 (Af p ) is a function of n p through its dependence on (f>(Af p ), 



M 2 o 2 M 2 
iPHK) = iPLn(n p ) = — - ^-^Af p {n p ) . 



(C5) 



Now we trace back the evolution to earlier times before the end of transition, Af > Af p . For this stage, instead of 
Af, it is more convenient to use n which is the number of e-folds from the critical point counted forward in time, ie, 
n = n p + Af p — N. Then tp 2 (n) is given by 

^ 2 (n) = exp2[/(n) - f {n p )]i> 2 min {n p ) , (C6) 

where f(n) a function that describes the tachyonic growth of the waterfall field, as given by Eq. (6.5). During this 
era, as for 0, we have 

0(n) - c e-™ . (C7) 
Here we note that n depends on n p and Af in a rather non-trivial way, 

n(n p ,Af) =n p +Af p (n p )-Af. (C8) 



In particular, we have 



where we have used the fact that 



^ = 1 + ^ = 1-2^, 
an p an p 



dAf, 



dn 



= —2rAf, 



V ! 



(C9) 



(CIO) 



which follows from Eq. (C4). 

Keeping in mind the above dependence of n on n p and A/", let us take the variation of (C6) and (C7). We obtain 



5ip 2 {n) = -2i}j 2 (n)f{n)5Af 
+^ 2 {n) (2 



f'^-nnp) 



d^ 2 min {n p ) 
dn„ 



8n p . 



dn 



i) = —r(f)(n)5n = rcf)(n)5Af — r<j>{n)— — Sn p . 



(Cll) 



(C12) 



By using Eq. (C5), one finds 



dn„ 



9M 2 r ^ 
J 2 J ^rAfp<2r«l 



where we have used Eq. (C2) and the fact that <j> p < <fi c in the last step. 

To find the variations of the number of e-folds, one solves Eqs. (Cll) and (C12) for 5Af to obtain, 



^ 5^p 2 (n) dn 1 5(j){n) 
tp 2 (n) dn p -2f'(n p ) r<j){n) 



2f(n) dn 



-2f'(n p ) dn r _ 



(C13) 



(C14) 



in which we have ignored the term \d\nip 2 (n p ) /dn p \ ~ r <C 1. Note that as a result of the vacuum domination 
condition, we also have \dAf p /dn p \ = 2rAf v <C 1. Then the above reduces to 



6Af= -■ 



1 <5V 2 M 



f(n) 



f'(n p ) 



5<p(n) 
r<j)(n) 



(C15) 



2f'{n p ) yj 2 (n) 

As clear from this result, the curvature perturbation spectrum is smoothly taken over from 8(j) to Sip 2 as n increases 
from to rip. 
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